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Abstract 

Based on representations of the symmetric group Sn, explicit and exact 
Schrodinger equation is derived for U = 00 Hubbard model in any dimensions 
with arbitrary number of holes, which clearly shows that during the movement 
of holes the spin background of electrons plays an important role. Starting 
from it, at T = we have analyzed the behaviour of the system depending on 
the dimensionality and number of holes. Based on the presented formalism 
thermodynamic quantities have also been expressed using a loop summation 
technique in which the partition function is given in terms of characters of 
Sn. In case of the studied finite systems, the loop summation have been taken 
into account exactly up to the 14th order in reciprocal temperature and the 
results were corrected in higher order based on Monte Carlo simulations. The 
obtained results suggest that the presented formalism increase the efficiency 
of the Monte Carlo simulations as well, because the spin part contribution of 
the background is automatically taken into account by the characters of Sjv- 
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I. INTRODUCTION 



The physics of holes moving in a background build up from a great number of electrons 
has received much attention in the last decade given by the connection of this problem 
with main questions arising from the study of strongly correlated systems (Chernyshev and 
Leung 1999). Developments in several fields underline this aspect: behaviour in antiferro- 
magnetic background (Anderson 1997), metal-insulator transition (Trugman 1990 a), doping 
effects (Loius et al. 1999), ferromagnetism (Kollar et al. 1996), high T c superconductivity 
(Trugman 1998), colossal magnetoresistance (Ishihara et al. 1997), spiral states (Watanabe 
and Miyashita 1999), generalized-statistics (Long and Zotos 1992), pairing mechanism and 
bound-states (Trugman 1990 b), polaronic effects (Louis et al. 1993), are main examples to 
be mentioned. From the theoretical side the Hubbard type models are especially involved 
in the theoretical description starting from infinite on-site repulsion (Nagaoka 1966), finite 
Hubbard interaction (Sutto 1991), and continuing with the inclusion of spin-independent 
interactions (Tasaki 1989), extensively extended versions (Kollar et al. 1996), degener- 
ated orbitals (Shen and Wang 1999) and limiting cases like the t — J model in one-band 
(Chernyshev and Leung 1999), or multi-band versions (Shen and Wang 1999). The theo- 
retical descriptions started with the 1-hole problem discussed by Nagaoka (Nagaoka 1966), 
which has been continued with two-holes studies (Kuz'min 1993), and many-hole analyses 
in various configurations and circumstances (Mattis 1986, Doucot and Wen 1989, Barbieri 
et.al. 1990, Toth 1991, Tian 1991). Despite the huge intellectual effort spent in the field 
exact results related to the mentioned problem are extremely rare (Mattis 1986) and con- 
cern only the 1-hole Nagaoka-state (Nagaoka 1966) and its possible extensions (Tasaki 1989, 
Siitto 1991, Toth 1991). However, in the above enumerated concrete applications mostly 
multi-hole descriptions are needed mainly in strong-coupling limit. These situations have 
been treated in the literature based on various approximations (acceptable in restricted 
conditions) like complete separation of charge and spin degrees of freedom (Kuz'min 1993), 
extended Gutzwiller wave-functions (Gebhard and Zotos 1991), canonical transformations 
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(Chernyshev et al. |cond-mat 98060 18|) , generalization of ID results to higher dimensions 



(Doucot and Wen 1989), restrictions related to the number of holes (Barbieri et.al. 1990), 
unrestricted Hartree-Fock (Loius et al. 1999), etc., or numerical procedures like exact diago- 
nalization on small clusters (Takahashi 1982), and diagonalization within a retained portion 
of the Hilbert space (Trugman 1998). We witnessed also studies on specific problems like 
ferromagnetism (Vollhardt et al. 1998). 

On the background of the mentioned large spectrum of description attempts, in this paper 
we are presenting an unifying aspect of the hole-motion in the simple one-band Hubbard 
model in the limit of infinitely large on-site repulsion U. For this reason first of all we 
have deduced and presented an exact and explicit Schrodinger equation for an arbitrary 
number of holes in the analyzed model. This equation clearly underlines the role played by 
the charge and spin degrees of freedom during the movement of the holes in the U = oo 
Hubbard model, being mathematically based on the representations of the symmetric group 
S*7v of degree N, where N represents the number of electrons within the system. Advantages 
of the presented description are underlined by the observations that: a) the emerging matrix 
elements can be find (or deduced in case of larger systems) from standard textbooks of the 
group theory, and b) the analysis reduced the concrete equation to be solved in subspaces 
fixed by a given, but arbitrary total system spin S of the initial large Hilbert space, allowing 
comfortable numerical applications. 

After this step, and concrete applications of the deduced equation related to the T = 
spectrum and physical properties of small samples, we turn on the study of T ^ properties 
by showing how the presented formalism helps the deduction of the partition function Z 
based on the characters of the representations of the symmetric group (again accessible 
from standard text-books in group theory). The emerging coefficients in Z can be reduced 
to coefficients easily obtainable from numerical methods, being related to the number of 
different loops with fixed length that can be drawn in the studied system connecting its 
lattice sites. The procedure is able even to increase the efficiency of the usual Monte Carlo 
simulation methods. Starting from the deduced formulas, we are presenting some examples 



of the T^O behaviour of the two-hole case. 

Our results clearly underline that in the infinitely repulsive Hubbard model the movement 
of holes is intimately influenced by the spin-background of the electrons. As a consequence, 
in describing this process in the general case, the charge and spin degrees of freedom cannot 
be entirely separated. 

The remaining part of the paper is organized as follows: Sec. II. describes the studied 
model and introduces the used notations, Sec. III. presents the basis of the Hilbert space 
used during the paper, Sec. IV. describes the action of the Hamiltonian on the basis vectors, 
Sec. V. presents the deduced Schrodinger equation together with applications, Sec. VI. 
concentrates on T 7^ properties and emphasises the advantages of the presented formalism 
in deducing thermodynamic quantities like partition function, thermodynamic potential, 
specific heat, expectation value of the total spin, Sec. VII. presents a summary of the paper, 
and an Appendix containing the mathematical details closes the presentation. 

II. THE MODEL AND NOTATIONS USED 

Our starting Hamiltonian (H) describes the one-band Hubbard model in U = 00 limit. 
Given by the absence of doubly occupied sites in the lattice, H becomes 

£«, = £ t tJ Hg = J2t tJ P ( £ cl c ha )P, (1) 

where t^- = — t are hopping matrix elements for nearest neighbour sites (otherwise = 0). 
The operators cJ jCr ( Q jCr ) creates (annihilates) an electron with spin a at lattice site i, and 
hi — J2a represents particle number operator. The double occupancy is projected out by 
P = (1 — nil )■ F° r convenience, we denote the lattice by A, the number of the lattice 
sites by N\, the number of electrons by N, so the number of holes becomes N h = — N. 
Being also interested in magnetic properties, we present below the notations and definitions 
related to spin characteristics. The total spin operator is taken as 

s = ^EE C w V (2) 

ieA cr.cr' 

4 



where r = (t x ,t v ,t z ) are the Pauli matrices. Furthermore we have S ,± = S x ± i S y and 
£ 2 = ( S x ) 2 + ( S y ) 2 + ( ) 2 . We denote by S z the eigenvalue of S z and by 5 ( S + 1 ) 
the eigenvalue of S 2 . The quantity S 1 represents the total spin of the state described. 

III. THE HILBERT SPACE VECTOR-BASIS USED 

Now we prepare a basis for a sector of the Hilbert space of the system with a fixed Nh 
number of holes and a given S z projection of the total spin. We start from Hilbert space 
vectors characterized by the position of holes h\ , ... , hN h , and the spin configuration of the 
electrons situated in single occupied states {0i}j 6 A\{/u,...,/ijv }> defined as 

i^w,w)= n c«i°>» ( 3 ) 

ieA\{hi,...,h Nfl } 

where the product is taken over all the occupied sites in A, and |0) represents the bare vacuum 
with no fermions present. For mathematical convenience, we introduce the convention that 
the elements of A are building up an ordered set (i.e. the lattice sites are numbered) and the 
product in Eq.(^[) follows this order. As a consequence, the product of creation operators 
from Eq.(|3p is written in increasing order with respect to lattice indices. Let us introduce 
now the function TZ : { 1 , . . . , N& } — > { 1 , . . . , N\ } that performs a regroupation of lattice site 
positions, maintaining however their initial numbering. By definition, 1Z gives the indices 
of occupied lattice sites in increasing order for the first N numbers in the {1, ...,Na} set, 
and introduces the indices of the empty sites in arbitrary order for the last Nh positions of 
the mentioned set (we underline, N\ = N + Nh). The numbering of hole positions obtained 
after the action of 71 being arbitrary, 1Z itself is not uniquely defined. Using it, \ip{h},{a}) 
becomes 

N 

\m,{«}) = n 4(*w w i o ) . (4) 

i=l 

In order to describe the spin configuration with a fixed S z projection of the total spin, 
we observe that S z can be characterized by permutations V G of the symmetric group of 
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degree N. Indeed, the set { Oi }ieA\{hi,...,hN } represents in fact a permutation of iVj- symbols 
| and = — 2 S z symbols I (see Appendix). Using "P, the spin configuration can be 
denoted as (ai)j e A.\{Ai,...Aiv } = ( °"( * ) ) )i=i> where the function a defined by cr( i ) := | 
if i = 1 , ... , iV| and u(i) :— J, if % = iVj + 1 , . . . , iV. With the introduced notations 
the basis vectors of the Hilbert space with fixed number of holes and fixed total spin 
projection S z can be written as 

N 

I K , V > = II 4 (i) ,^») I ) • (5) 
1=1 

The subspaces generated by vectors { | 1Z , V ) }pes N with different 7Z are different but 
isomorphic parts of the whole Hilbert space of the system. Let us denote by 7Yf z the Hilbert 
space which is isomorphic with all subspaces generated by {\7Z , V)}v&s N with arbitrary but 
fixed 7Z and S z (see Appendix). Similarly 7i c denotes the Hilbert space which is isomorphic 
with all spaces generated by vectors 1 1Z , V) with arbitrary but fixed V and S z . We mention 
that independently on S z , 7i c represents in fact a Hilbert space of a system of hard core 
spinless particles. With this notation the whole Hilbert space of the system TC is a direct 
sum over S z sectors Ti = ( 7i c <g> 7Yf z ). We underline, that introducing the notation 
7i s = ©5*Wf , we have only H = H c <g> W s . 

In the remaining part of the paper we denote an arbitrary element of Hf z , 7i s and 7i c 
by (/?f z , and <^ c , respectively. Furthermore the basis vectors of 7i c will be written as 1 1Z), 
and the basis vectors of 7i s s z as \V). With these notations the basis vectors that we have 
are \7Z , V ) = \7Z) <g> \V ). The states | 7Z ) describe in fact the charge distribution 
and | V ) the spin configuration. We have to mention that Eq.([|) has to be used with care 
because it contains redundant information. Every 7Z G Sn a and V G Sn uniquely define 
a basis vector \ 7Z) ® \V), but several different pairs of permutations define the same 
basis vector. Because of this reason indexing these basis vectors by integers, in the case of 
the basis set { | IZi } } C 7i c the index i goes up to {^n and not to N !. Similarly, the 
basis set { | V% ) } C 7i s s z has Lyv 2 J + s z ) different elements. The reason why we express the 
basis vectors in term of permutations is detailed in the Appendix: the Hilbert space of the 
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spin configurations for fixed N and S z is isomorphic with a proper right ideal of the group 
algebra C[ } of the symmetric group Sn- 

We emphasize that the Hamiltonian presented in Eq.(|I|) cannot be written in the product 
form 7^ H c ® H s and the effect of the Hamiltonian on 7if z depends on the hole 
configuration (charge distribution) 1 TZ). In fact we have HQ = HQ ® HQ (TZ), therefore 
the eigenstates of H in general case cannot be written in the form | ip c ) <g> | <p s )• As a 
consequence, there is no apriori charge spin separation taken into account at the level of the 
description. However, contrary to H, the effect of spin operators does not depend on \7Z), 
i.e. on Tl c . This is why the decomposition of the Hilbert space into sectors characterized 
by a fixed total spin S can be done within Tt s = ©g* (®s=s* ^-T' S )- Denoting by (p sz ' s 
an element of 7i sz,s = T~L C ® Tt^ z,s ' , an arbitrary common eigenvector of S z and S 2 in the 
sector of fixed number of holes can be written as ip c <g> ' . 

IV. ACTING WITH THE HAMILTONIAN ON BASIS VECTORS. 

The contribution HQ of the Hamiltonian from Eq.([l]) has non- vanishing effect on a 
basis vector from Eq.(g|) only if i represents a hole position and j an electron position, i.e. 
1 < 7Z~ 1 ( i ) < N and N + 1 < 7Z~ 1 (j ) < N^. Furthermore, regarding the spin sum in 
Eq.(P, only the term c\ a , djy of HQ has non- vanishing effect in which a' = a(VTZ^ 1 (j ) ). 
Using the notation J = 7Z~ 1 (j ), we obtain 

N 

HQ ( I TZ ) <8> | V ) ) = %(j{v{])) CjMPO)) IT ^n(i)MV(i)) I ) = %a(V0)) A?>CP(D) %(i), CT (-p(i)) 

i=i 

■ ■ ■ Cn{3-i),a{v(3-i)) c j,a(vm c n{]+i)MV(5+i)) • • ■ c k(n)mv(n)) I ) 

= 4l(l),a(V(l)) ■ ■ ■ C n(j-l),a{V{3-l)) 4,a(P(j)) C ^(J+l),a(P(J+l)) • ■ ■ 4l(N),a(V(N)) I ) 
N 

1=1 

where V^> represents the transposition that interchanges the indices i and j. The result 
cannot be written in the form | TZ) ® | V ), because the permutation TZ does 
not meet the requirements regarding the site ordering presented after Eq.@. Indeed, the 
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creation operators are no more situated within the product in increasing order with respect 
to their lattice indices because ct ,.p,~^ is situated in a wrong place. To overcome this 
situation, let us introduce the cyclic permutation C( 1Z ; i , j ) G SV which rearranges 
the N pieces of operators of the last row of Eq.(|(J) in the required order by moving the 
iaCPQj)) °P era tor positioned at site j to the required position within the product. In the 
i < j (i > j) case, for this purpose C( 1Z ; i , j ) interchanges the 4o-(-p(j)) operator with 
previous (following) operators. The last operator with which we interchange, is situated in 
the position 7Z~ l (i' ). Here, in the i > j case i' is the least index among the indices greater 
than i for which 1Z~~ 1 ( i' ) < N, and in the i < j case i' is the greatest index among the 
indices smaller than i for which 1Z~ 1 ( i') < N. In other words i' labels the occupied lattice 
site whose index is closest to i among the indices situated between i and j. Now we have 

C(K;i,j) = C^- 1 ^- 1 ^')) . (7) 

It can be observed that we get the ordered sequence of the operators in Eq.(|6|) if we permute 
their indices by C _1 ( 1Z ; i , j). Given by this, also a sign-change ( — 1 y c ( n ' l M is obtained, 
where | C \ has the meaning of the parity of the permutation. Based on the presented 
properties we obtain 

N 

ffij M 7? \ 6d \ V\ ) — ( — 1 )\ c ( n ^^)\ TT r 1 I n \ 

JJ ooVl /v / w I ' I) \ L ) 11 t -7'(i.j)7eC- 1 (7Jii ) j)(/) ! o-CPC- 1 (K;i,j)(0) 1 ' 

1=1 

= \ V^ ) TZC-\TZ; i,j)) ® \(-l)W R -M\VC-\K\ (8) 



The difference between the hole configuration described by the permutation 1Z and 
j>(ij) C~ 1 ( 1Z ; i , j ) lies in the position of one hole, which has been moved from site 
i to site j. This shows that the holes are moved by the Hamiltonian as if they were hard 
core bosons, while the spin configuration is changing depending on the actual hole config- 
uration. The Hilbert space generated by the vectors | 1Z ) is actually the Hilbert space of 
Nh hard core bosons (| JZ ) describes the state in which TV/, particles are situated in the 
1Z( N + 1) , . . . , 7Z(N^ ) lattice sites positions). Therefore, we write the effect of the Hamil- 
tonian as -f/fcoo 1^)5 taking into account supplementary the effect on the spin configuration. 
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According to the Appendix, the effect on the spin configuration can be described with a 
linear operator T SZ [C( 71 ; i , j ) ] where T s * is a representation of Sn uniquely determined 
by S z . Based on these considerations, the effect of the Hamiltonian on the basis defined by 
Eq.(0) can be given as 

Hoo{\K) ® | V )) = £ty#gjft) ®T s *iC(1l;i,j)]\V). (9) 

Starting from Eq.(JD we analyze below the Schrodinger equation describing the behavior of 
the system. 



V. THE SCHRODINGER EQUATION 

Due to the fact that the system is invariant under the SU(2) global rotations, in the spin 
space the Hamiltonian and S 2 are simultaneously diagonalizable. Therefore the spin eigen- 
subspaces are invariant under the effect of the Hamiltonian. The effect of the Hamiltonian 
on Ti,g Z ' s can be described by a T$ representation of the symmetric group independently of 
S z (see Appendix). Let us denote by (Ts[V) )„ m the matrix elements of the operator Ts[V] 

in an orthogonal and normalized basis of 7i^ z,s with arbitrary S z (denoted by { }), and 

i ■ ~i- ( Na ) 

by ( .f^oo )ki = ( T^h I -f^oo I T^i ) where IZk , Tli E { \ IZi ) }}J{ are basis vectors indexed 

by integers in T~C C . Now the Schrodinger equation becomes 

a lrn Yl f ij ( ^oo )kl ( T s [ C( Tli ; % , j ) ] ) nm = E a kn (10) 

l,m i,j 

Eq.(|10) gives the energy eigenvalues E connected to the eigenstates with given total spin S. 
The akn coefficients in Eq.(^) represent the coordinates of the energy eigenstate E in the 
basis 1 1Z ) <S> \ b s ) given as 

foM = £ a *J^*> ® IK) (ii) 

fen 

We have to mention that essentially the same wave vectors as presented in Eq.(|TT|) has been 
used by Long and Zotos (Long and Zotos 1993) in the study of the stability of the Nagaoka 
state, without however to write an explicit Schrodinger equation for it. 



It can be seen from Eq.([TT|) that the matrix elements of the Hamiltonian between spin 
states depend also on the hole configuration TZ, i.e. charge distribution, so dynamically the 
spin degrees of freedoms are not separated from the charge degrees of freedom. At the level 
of the basis states however, we have the possibility to treat them separately and this allows 
us to write the Schrodinger equation in the form of Eq.(|Tt]). We also can use Eq.(|TD|) for 
exact numerical diagonalizations. The main advantage of Eq. fllPf) is that the matrix elements 
of the starting Hamiltonian have been reduced to a direct sum of matrices related to the 
invariant subspaces Hf sz . Because of this reason, in concrete applications we need to solve 
the eigenvalue problem of Eq . (|T0|) in a relatively small Hilbert space instead of 7if z which 
is much larger. For concrete application the following indications are needed in order to use 
Eq.(|10]): The matrix elements of the hard core boson Hamiltonian can be easily calculated 
and for small number of electrons the matrix elements of the representation Tg is given in 
mathematical textbooks (see for example Hamermesh 1962, Table 7.3). For greater N values 
their clear calculation procedure is also available from same sources (see e.g. Hamermesh 
1962, Section 7.7). 

Eq.(|Tt]|) shows explicitly the role played by the total spin S in the Schrodinger equa- 
tion. For example in the case of S = S max , the representation Ts max is the alternating 
representation. In this case the dimension of 'j-i^^^-^ is 1, and the matrix representing 
every even permutation is 1 and every odd permutation is — 1. Because of this reason hard 
core fermions with maximal spin do not differ from spinless fermions (see also Long and 
Zotos 1993). The on-site interaction has no effect on spinless fermions therefore there is no 
differention between hard core and free particles in this case. 

Another simple example is the ID case with open boundary conditions. Numbering the 
sites of the chain in order one after another, the form of every nearest neighbour pair is 
simply ( i , i + 1 ) in this case, and we have C(TZ ; i , i + 1) — 1 G SV whose matrix is 
the unity matrix for arbitrary representation. Therefore there is no difference between hard 
core bosons and fermions, so the energy spectrum is independent of the spin of particles, 
only every energy eigenvalue is ( 2 S + 1) N fold degenerate, S being the spin of particles. 
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If D > 1 and S z 7^ .Smaa; than the effect of the Hamiltonian on the spin background 
is highly nontrivial, and influences the mobility and interaction of holes. A such type of 
situation was analyzed by numerical technique by Trugman (Trugman 1998). He showed 
that on Neel type antiferromagnetic background one hole is mobile but two holes are less 
mobile than it was previously considered. This is due in fact in the language presented here 
by the permutation possibilities of spin positions with mobile hole positions. Trugman's 
results underlined as well that the effect of the Hamiltonian is clearly seen also on the spin 
background, which cannot be treated entirely separated from the movement of the holes 
within the system. 

We note that the Schrodinger equation Eq. (^) essentially differs from Kuzmin's equation 
(Kuz'min 1993) for the purpose to solve the two hole problem on singlet background. In the 
mentioned reference the Schrodinger equation has been solved supposing implicitly that the 
Hamiltonian has no effect on the spin degrees of freedom i.e. spin background. From Eq.flTOl) 
it can be seen that this is correct in the case of a hard core boson system only. This explains 
why Kuz'min obtains the result that the Nagaoka state (the fully polarized ferromagnetic 
state) is degenerated with the singlet state even in the one hole case, in contradiction with the 
Nagaoka's original result (Nagaoka 1966). In one hole case the spectrum given by Kuz'min 
is really related to a hard core spinless boson system spectrum, which is equivalent with 
the spectrum of a spinless fermion system. This last one, on its turn, is equivalent with 
the S max spectrum of our original electron system as presented above. The first equivalence 
is due to the fact that for hard core spinless particles (both bosons and fermions) the N 
hole problem is equivalent with the N particle problem, and in the one particle case the 
statistics is unimportant. For more than one hole it is not the same whether we deal with 
bosons or fermions even for hard core interaction. The ground state energy of the bosons is 
always smaller (Long and Zotos 1993). In order to exemplify, we mention that the equality 
presented by Kuz'min 



o = E 



1 



(12) 



k 



E — ( €k + tp-k ) 
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gives the spectrum of a spinless hard core boson system consisting of two particles (or 
equivalently two holes). Here is the one particle dispersion relation and P is the total 
momentum of the system. For the ground state energy of the simplest non-trivial two 
dimensional case, the 2x3 lattice, Eq.(|l2"D gives Eq = —6.6468 while from Eq.(|H]) one 
find Eq = —6 which connected to the S = S max = 2 value of the total spin, and 
££=° = -5.0212 is the lowest energy value on singlet background. Furthermore we note 
that in the two particle case the spectrum of hard core spinless boson system is identical 
with the spin 1/2 fermion system's singlet spectrum, therefore it can be computed using 
Eq.flljD. However, as presented here, in the S ^ S max case the behaviour of the system 
described by Eq.(|l|) cannot be described by hard core particles and independent spin degrees 
of freedom. This is the reason why the iV particle and the N hole problems are not generally 
equivalent, and the behaviour of few holes is more complicated than the behaviour of few 
particles. 

VI. THERMODYNAMIC QUANTITIES 

An asset of the representation presented in Eq.@ is that the trace of ( ) l over the 
spin degrees of freedom can be computed with this relation exactly. This allows us to 
calculate expectation values at T ^ of different thermodynamic quantities. The procedure 
is presented as follows. 

A. The partition function and total spin 

If we are interested in the magnetic properties of the model Eq.(p]), we may express for 
example the expectation value of the square of the total spin. We have 

(S 2 ) = WS(S + 1) £ Tr^e-#4E3(^) 2 Tr w , e -#, (13) 

Z S -S<S*<S Z S z 

where the partition function is given by 
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E E Tt h^s e "t = E Tr W 5^ e "t . (14) 

S -S<S Z <S S z 



In Eqs.( |l3"| , |l4"D we used the fact that Tr n s z ,s e kr does not depend on S z due to the SU(2) 
symmetry, therefore Y^s= a s z Tr w s z ,se ^ = Tr w s* e~~^ for arbitrary S*. As a consequence, 



the subspace under the Tr in the last relations from Eqs. (^L3| , p~4|) means that the trace should 
be taken only over the Ti sz subspace. 

For finite lattices the Tr operation means a finite sum which can be interchanged with 
the sum arising from the series expansion of the exponential function. We need in this case 
for calculations Tt h s* ( Hoo ) l . This trace can be computed starting from Eqs.(§, [|). We 
have 



H3l X 



Tr^iH^y = EE E ••• E E •••E^ 1 

\v) \n) 'ui k,ii 

( ft | | TZw ) ... (Ki\Ht£\K) x 

(P|T 5z [C(^_ i; i,, j,)] ...T sz [C(n 1 ;z 2 ,j 2 )}T sz [C(n; ii,h)]\V). (15) 

Let us denote for 7 (/) = 7 (0) by 7 = ( 7 (0), 7 (1),..., 7(^-1), 7(0) = {TZ , TZ X , . . . ,TZi- X ,TZ) 
a sequence of TZ permutations for which T\ i=1 (7(01 -^fc,oo 1 1( i — 1 ) ) 7^ 0. We call L( 7 ) := / 
the length of 7. The set of all the 7's with length I will be denoted by Qi. Let consider now 
the lattice consisting of every different hole configuration TZ. We recall that TZ can be 
described by the ordered position of holes: (hi < ... < h Nh ). We consider two lattice points 
TZ and TZ' nearest neighbours if they differ by only one hole position, and these different 
hole positions are nearest neighbours in the original lattice A. The lattice that we obtain 
in a such a way is a part of the D Nh dimensional hypercubic lattice, and 7 is a sequence 
of nearest neighbour lattice points in A^. Because of this reason, we call 7 a ,,loop". In 
Nh = 1 case 7 is a loop in the original lattice A. 

In the case of Ya,j Uj ( TV I | TZ ) ^ there is precisely one nearest neighbour pair 
( i , j ) for which ( TZ' \ \ TZ) ^ 0. Therefore the loop 7 uniquely determines a sequence 
of pairs of nearest neighbour indices ( i\ , ji ) , ... , (ii , ji) which gives the only one nonzero 
contribution of the sums S^jj • • • Eijjj i n Eq.(|l^). Denoting the product of permutations 
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C{ 1Z ; z , j ) obtained from this nonzero term by P 7 , we have 

Tr^^y = £ EW^t^lW = E x S2 (^)- (16) 

As can be seen from Eq.(|TJ|) we obtain a sum of characters of T sz [ V y ]. Every character is 
constant on an arbitrary conjugate class C. Therefore, the above character sum has N^ C >(1) 
identical members, where ( I ) represents the number of the loops with length I for which 
Vry G C. Using the notation f3 for t-times the reciprocal temperatures and the results for 
the sum of the characters x S * from Eqs. ( p6] - |27|) from the Appendix, we obtain 

(s 2 ) = MnE a^>(0(-i) |c| (2^) (fya) , (17) 

4 z 1=0 11 ccs N v / V i=i / 

where, for the partition function we have 

z =£7! E iv (c) (/)(-i) |c| (tiZ**) ■ (18) 

1=0 1 ■ ccs N v j 

Here | C \ is the parity of permutations from the conjugate class C and ( C\ , C 2 , . . . , Cjv ) 
describes their cycle structure, i.e. these permutations contains Cj cycles with length i. 

Concerning the technical aspects in using Eq.(^) or Eq.([i7|) in concrete applications 
we mention that in order to deduce the coefficients N^ c \l) one should follow the following 
procedure. We need to deduce a concrete J2ccs N contribution at a fixed / value. For this: 
1.) Start from a fixed hole configuration 1Z\, and denote with different numbers every site 
occupied by electrons. 2.) A single step for a given hole it means that we have to interchange 
the given hole with a nearest neighbour number. 3.) Take I steps with the holes in such a 
way, that finally get back the original hole configuration (this is a loop with length I). In this 
process an arbitrary hole can be moved in every step, the order of the steps being relevant. 
However, by interchanging holes between them, we do not obtain new loops. 4.) Determine 
the cycle structure C of the permutation of the numbers for the obtained loop (i.e. obtain 
the numbers Cj). As a result of this analyses, we have find in this first step m = 1 one 
loop of length / and a given cycle structure {Cj} m= i. 5.) Go back to the step 3). and find 
a different loop with the same length I based on the same starting hole configuration TZ\. 
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After this step m = 2, we have another loop with length I, and another cycle structure 
{Ci} m= 2- 6) If you have taken into account every possible different loop starting from the 
hole configuration 1Z\, do the same procedure with every possible different starting hole 
configuration TZj. 7) All this being done, the summation J2ccs N a t a fixed I means in fact a 
summation over all possible cycle structures obtained above. For example, given a concrete 
cycle structure {Ci] m = {C\ = 01,62 = 02, ...,C r = a r ) we may express with it J2iCi or 
J2i i 2 Ci, the quantity N(°'(l) having the meaning of how many times this concrete {Cj} m 
cycle structure has been obtained in the procedure presented above. The number | C | gives 
the parity of all permutations which give the same cycle structure {Ci} m (as presented in the 
Appendix, all different permutations with the same cycle structure have the same parity). 
8) As presented in the step 7), we have to sum over all possible cycle structures obtained in 
points l)-6). 

Using this procedure, with a computer, a partition function or a spin-square T ^ 
expectation value can be deduced for a system with arbitrary number of holes. Evidently, 
the calculation time increases with the size of the systems, or number of holes. 



B. The free energy and specific heat 



^From the partition function Z given in Eq . (|18|) we can compute the free energy F = 
— T In Z and the specific heat c = T 9 ^g^z Z ^ ■ 

From Eq. ( |i~7l) we correctly obtain paramagnetic behaviour at high temperatures in ther- 
modynamic limit. The first term of the /-sum becomes dominant if T increases, and 
N( c '( ) = 1 if C = 1 otherwise N^ c \ ) = 0. Therefore we have in high temper- 
atures limit ( S 2 ) = 3 / 4, so the value of the total spin per particle is proportional with 
1 /y/N. 

In order to analyze the temperature dependence of ( S 2 ) we need the values of ( I ) 
as functions of I. Comparing Eq.flTB]) with the formula for the partition function Z derived 



in first quantized formalism we obtain 
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— J : — - = i-^-i (19) 

for arbitrary conjugate class C which can occur in Eqs.( |iT| , |TE| ), where Eq is the ground state 
energy of the N particle hard core boson system. 



C. The Nagaoka state 

The knowledge of the asymptotic behaviour given by the Eq. ([[£]) is enough to analyze 
the T dependence of the Nagaoka ferromagnetism (i.e. = 1 case). In this situation Qi 
contains the I length loops of the lattice A. For one hole, 1Z is uniquely determined by the 
position of the hole. The positions of the hole is described by r y{ ) at the starting point of 
the loop, and 7( I ) = V n J ' . . . V n J1 7( ) V^ 1 at the ending point of the loop are the same. 
It can be seen that the parity of | V 1 \ = ( — 1 ) l is always even when we use open boundary 
conditions, because there is no loop with odd length in this case. It means that the dynamical 
evolution controlled by the Hamiltonian Eq.(|l|) is not able to permute the particles by odd 
permutations, given by the presence of the hard core potential acting between them. As a 
consequence, the fermionic character of the particles has absolutely no effect in this case. 
This is fully consistent with our previous statement (presented before Eq.(|l^) ) that the 
S = S max spectrum of our original fermion system is equivalent with the spectrum of a 
hard core spinless boson system. Moreover, if bosons ,, existed" with half spin, then there 
would not be differences between the system consisting of these particles with hard core 
potential on a square lattice and our original fermionic system in the one hole case, because 
the permutations C( TZi ;i , j ) occur in the Schrodinger equation Eq.flTUD are always even, 
and for even permutations there are no difference between the representations T$ connected 
to bosons and fermions. As we mentioned above, the ground state wave function of a boson 
system is always symmetric, therefore the spin wave function is also symmetric, i.e. the 
ground state is ferromagnetic. This extreme "symmetry" is an interesting explanation of 
Nagaoka's theorem. Our formalism certainly gives back this result. It is obtained by taking 
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the T — > limit in Eg . (|17D , considering the asymptotic behaviour of the coefficients N^°\l) 
given by Eg . (|19D for the conjugate classes with even parity and N^(l ) = for odd parity, 
furthermore using the result Eg . fl25|) from the Appendix. When we use periodic boundary 
conditions we still have \V y \ = ( — 1 which, however, can be odd, but only in the case 
when the linear length of the lattice is odd at least in one direction, and 7 goes through the 
boundary in this direction. In t < case (3 l ( — 1 is always positive, the sign coming from 
the odd parity of the permutation is compensated by the sign obtained during the movement 
around odd steps, therefore the above argumentation remains valid. The key feature now is 
the fact, that the parity of V 1 is in one to one correspondence with the parity of L( 7 ). 

D. Cases with more than one hole 

The existence of more than one hole permits arbitrary (even and odd) permutations 
of the particles under the dynamical evolution (independently of the length of the loop), 
therefore the fermionic feature of them play an important role in building up the energy 
spectrum of the system. Using the asymptotic behaviour of N^ c \ I ) from Eg. (pi]) and the 
results of the Appendix presented after Eg. (|25|) , the ( S 2 ) value remains undetermined (i.e. 
0/0), which is due in our interpretation to the fact that Eq is not a possible energy for the 
fermionic system any longer. The special similarity between bosons and fermions does not 
remain valid, therefore we can not draw conclusion for the magnetic property of the system 
based on Nagaoka's theorem. There is a gualitative difference between the one and two 
hole case. It is similar with the case when in a special point the Hamiltonian has an extra 
symmetry, but the features of the model could be very different even close to this point 
due to the symmetry breaking. If bosons would have half spin, a system consisting of such 
hard core bosons on a lattice would be ferromagnetic even for arbitrary hole concentration. 
This statement remains valid for real bosons with integer (non zero) spin as well and the 
expectational value of the total spin goes to a macroscopical (but not to the maximum) 
value when the temperature goes to zero, in this case. 
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In order to obtain information about the two (or more) hole case, the study of Eq.(pj]) is 
needed. To deduce the ( I ) numbers one can use Monte Carlo methods. For example, 
by random sampling we can determine the percentage of the contribution of loops with a 
given length in V y G C . However, we have to mention, that the trace over the spin degrees 
of freedom has been taken already, and because of this reason, the presented procedure 
increases the efficiency of the usual world line algorithm. In this case, we need only a 
sampling from the world lines of the holes, which are 2 N times less than all the different 
world lines of the electrons. To see the trend of the magnetic behaviour in the T dependence, 
we have determined the ( / ) coefficients up to the 14-th order exactly, and further, from 
15-th up to the 50-th order by Monte Carlo in the case of a 5 x 5 lattice with two holes. 
The square of the total spin per particle and the specific heat deduced from these results are 
presented in Figl.a. and b. As the temperature decreases the error of our results increases, 
but it is clear that the total spin increases, which underline a tendency to ferromagnetism 
at T — ► within the accuracy of our calculation. 



VII. SUMMARY 

Summarizing, we have analyzed a simple one-band Hubbard model in U = oo limit 
in any dimensions. We have treated this problem for the case in which an arbitrary but 
fixed hole concentration is present within the system. Concentrating on the effect of the 
Hamiltonian on the charge and spin degrees of freedom and based on the representations of 
the symmetric group, we presented an exact Schrodinger equation in a form which explicitly 
shows the effect of the spin background in the movement of the holes. Explications related 
to the application of the presented equation in concrete situations were given in detail. 
Based on the deduced results we have showed that the effect of the spin background on 
the movement of holes is unimportant a) when the system is one dimensional with open 
boundary condition and b) in any dimensions when the spin of the background is maximum. 
In the first case the holes move as hard core spinless bosons, in the second case as hard core 
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spinless fermions. When we have D > 1 and S ^ S max the effect of the spin background 
seen in the motion of the holes and the energy spectrum is highly non-trivial. The Nagaoka 
ferromagnetism (A^ = 1 case) in our formalism is due to the fact that in the one hole case 
and U — oo, the rearrangement of particles with odd parity permutations in the dynamical 
evolution of the system is impossible. Therefore, the fermionic feature is not playing any 
role, and the system is equivalent with a system of hypothetical hard core half-spin boson 
system whose ground state is trivially ferromagnetic. In the presence of more than one hole 
this is no more true and the complete separation of charge and spin degrees of freedom is 
impossible. 

On the other hand, based on our representation of the Hamiltonian, the trace over the 
spin degrees of freedom can be computed exactly. The presented procedure allows us to 
express the partition function, free energy, specific heat and the expectational value of the 
square of the total spin. The calculations can be given based on the counting of loop 
contributions, whose steps were described in extreme details for potential applications. In 
principle the method allows to calculate T ^ expectation values for an arbitrary number 
of holes, based on clear formal steps handable by computer. In the frame of this formalism 
the efficiency of the usual Monte Carlo methods can also be increased. 

APPENDIX 

In this Appendix we are presenting mathematical details related to deduced formulas 
presented in the paper. For the beginning we are presenting an overview of notations and 
definitions used, then rather technical proofs follow. For more mathematical details con- 
nected to the deduction procedure we refer to Hamermesh 1962. 

We treat a permutation of degree A as a bijective function V : {1, . . . , A} — > {1, . . . , A}. 
The product of two permutations is defined by the standard composition of them as functions 
(from right to left): ( V Q )( % ) := V( Q(i) ). The symmetric group of degree A formed 
by all the permutations of degree A is denoted by S N . For M C { 1 , . . . , A } the group 
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Sm formed by all the bijections from M onto itself is a natural subgroup of Sn, a function 
of S M acts identically on the elements of { 1 , . . . , N } \ M. We denote the transposition 
which interchanges i and j by and C^ - ^) the cycle which moves i into j while pushing 

left (right) the elements between % and j with 1 if % < j (i > j). The group algebra C[Sn] 
is a complex Hilbert space generated by the elements of S N as an orthonormal basis. The 
associative but noncommutative product of two elements of the group algebra defined by 
the convolution: 

ab = ( ^ a v v) ( £ p Q Q) 

Ves N Qes N 

= E E a v (3 Q VQ= £ ( £ a Q P Q -i v )v. 

•PgSjv QgS n VeS N QeS N 

If a subalgebra S has the property that, for a in ab is also in S for any elements b of 
the whole algebra, then S is called a right ideal. If e is an idempotent element of the algebra 
(e e = e) then eC[S N ] = {ea|aG C[S N ]} is a right ideal, because of the associativity 
of the group algebra. It remains true if e is only essentially idempotent (e e = c e , c G C) 
because e' = e / c is idempotent. 

Now we verify that the Hilbert space generated by all the vectors \1Z , V ) with arbitrary 
but fixed 1Z is isomorphic with a proper right ideal of the group algebra C\ Sn ]■ (In 
the following iV and S z are fixed.) We recall that for fixed S z we denote by a the following 
function: a(i) :=|ifi E {l,...,iV T = iV/2 + ,S z } and <r(i) :=J,ifi G {iV T + l, . . . , N}. Let's 
denote the subgroup of Sn the elements of which permute the numbers of the set o" _1 ( j ) 
by X| = iSg-i (-]■). Similarly, the permutations of = S a -i(^ permute only the a _1 ( | ) 
numbers. Let K = x K^, the direct product of and K^. It can be seen that 

:= -±= E ^ (20) 

y I a | 

is essentially idempotent, (| K \ = ! N± ! is the order of the subgroup K) therefore 
c^s z ._ e sz C[Sn] is & right ideal. For different permutations V and Q the elements e sz V 
and e sz Q of the ideal ^s sz are the same if and only if the permutations are in the same right 
coset with respect to the subgroup K (i.e. V = V Q, V G K), because in the expression 
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J2veK VV Q = J2v"£K T"' Q — J^Q'eKQ Q' the sum is taken over the right coset K Q, 
and they form a disjoint cover of the group. 

Examining the basis vectors 1 TZ , V ) for a fixed TZ can be seen that the definition Eq. ([5]) 
gives the same vector for different permutations V and V V ( V G K ) from the same right 
coset, because a(V(i) ) = cr( V' V{ i ) ) , {i G { 1 , . . . , N } ). However, permutations 
V\ and V2 from different right cosets mix the numbers <J~ 1 { | ) and cr _1 ( | ) with each other 
differently, therefore they lead to different spin configurations. 

The above two paragraphs prove that the map | TZ , V) — > e V for fix TZ is well-defined 
and injective. Let's extend this map linearly over the whole space generated by the vectors 
I TZ ) ® I'P) with fix TZ. Certainly we get a surjective map onto Q sz , since the elements 
e s z p generate this ideal. Now we show that this linear isomorphism also preserves the inner 
product. In the proof we denote the inner product in the group algebra by ( . , . ). 

V Y I K I SieA' J I K I s 2 ei< / 

= T^l E E 2') ■ (21) 

If the right cosets KV and -ft" Q are different then they have no common elements, therefore 
( V , Q! ) = 0. (The elements T 3 G S'at are orthonormals by definition.) If the right cosets 
K V and K Q are the same then the sums have | K | different members which give 1 (in case 
V = Q'), the rest are zero, so the final result is 1. 

We showed that the Hilbert spaces generated by the vectors | TZ , V ) with arbitrary but 
fixed TZ are isomorphic with ( =s sz therefore also with each other. This space is denoted by 
Ti instead of ^s sz . The image of | TZ , V ) by this isomorphism is independent of TZ. This 
element of 7i, sz,s is denoted by | V). Based on the above defined Hilbert space isomorphism 
we identify Ti sz with Q sz and we distinct neither the notations of the inner product nor 
the operators act on them anymore. We can regard a linear operator defined on Q sz as an 
operator on Ti sz and vice versa. 

Similar technical and even more simplier work to prove that the space Ti c is isomorphic 
with the Hilbert space of the system of N (or equivalently Nh) hard core spinless particles, 
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and its basis characterized by the positions of these particles uniquely determines the basis 
vector | TZ ) . 
The formula 

T[Q]V := ( - l) |s| V Q" 1 (22) 

defines a representation of the symmetric group in the group algebra. Here T[ Q] is a linear 
operator associated the Q e Sjy by the representation T, and V is a vector as an element 
of C[ Sn ]. This representation is fully reducible. It is clear from the definition that every 
right ideal is an invariant subspace of this representation, thus T s * := T^s 2 is also a 
representation (with lower dimension) of the symmetric group in Q sz . Comparing the above 
defined representation and the effect of the Hamiltonian on a basis vector \1Z) ® | V ) as 
can be seen from Eq.(|||), it is clear that the effect on the spin configuration can be described 
by the operator T°* [ C( TZ ; i , j )], and we get Eq.(|9]). 

Now we compute the character of the representation T s * . By definition the character is 
the trace of the operator T s " [ V ] for a V G Sjy. 

X SZ {V) = £<2I^H^]|2> = -L E (e s * Q\T s °[V)\e s * Q) 

\Q) I A I Q£S N 

= r^r E (<Pq\{-i)W<Pqv- 1 ) = E E (shQi&Qv- 1 ) 

(_ 1)1^1 (-1) |P| 

= I f? 12 E E ^siQAGP- 1 = r&n E E ^q-pq- 1 ^ • (23) 

I I GeSV Si^aeJf I xv I QeS N SeK 

Every element V of the conjugate class C v which contains V appears as Q V QT 1 the same 

times. That's why there are | Sn \ / \ C v | permutations Q for which Q V QT 1 = V . (| X \ 

means the number of the elements of the set X.) If V' is in K then the second sum in the 

last row of Eq.(|2"B|) has one non-zero element, otherwise it does not. Therefore the character 

x sz (v) = (- 1)1*1 







c v nx\ 




K | 


\c v 





r> I r" ! r> I r" 



-1> -i j > -jv ' -JV 

Cj + c" = Cjf , . . . , c' N + = C 1 ^ 
c[ + 2c' 2 + . . . + Nc' N = N/2 + S z 
c'l + 24' + . . . + TVc^ = AT/2 - S z 
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n /r<r\ 

= (-i) |p| e n ; . ( 24 ) 

ci,...,c N i=i \ L i / 

Eti ic * = N / 2 + sz 

where the cycle structure of every elements of the conjugate class C v is described by the N 
pieces of numbers , , . . . , C N by the following way: these permutations contains 
cycles with length 1, cycles with length 2, etc. For the second equality see Hamermesh 
1962, the third can be obtained with simple computation. 

We need the sum of these characters on an arbitrary conjugate class C C where 
is the alternating group formed by all the even permutations of Sn- Based on the first line 
of the formula Eq.flSID we can see, that 



E \C\X S \C) = J" \A N nK*'\ = ^. (25) 

CcA N I K S Z I 2 

The similar sum over the conjugate class which has no common element with A^ gives the 
result — N \ / 2. The even permutations form an invariant subgroup in Sn therefore every 
conjugate class is in A^ or disjoint from it. Hence J2ccs N I C I X 5Z ( C ) = 0. 
Based on the last line of the formula Eq.(p4|) we can get: 



N/2 

E X S \C) = (-1)1°! 2^. (26) 

S z =-N/2 



Furthermore 

N/2 



N ' 2 v / 1 N \ 

E (S Z ) 2 X SZ (C) = (-1)^2^ . 

=-JV/2 V 4 i=l / 



(27) 

S*=-N/2 

The representation T sz is not irreducible if S z ^ S^ ax or S^ in . Certainly it can be 
verified in pure mathematical way, but physically it is due to the SU(2) rotational symmetry 
in the spin space. The operator S~ realizes a linear bijection between the spaces 7if z ' 5 and 
7ig Z ~ 1,s if S > S z . Moreover [ , S~ ] = 0, therefore the effect of the Hamiltonian 
is equivalent on the space TCg Z,s and 7if z_1 ' s . The effect of the Hamiltonian on the space 
Ji^max = ftSmax'Smax can described by the representation T s ™ ax =: Ts max , hence it 
is true for the space 7-t^™ ax ~ 1)Smax . However, the effect of the Hamiltonian on 7if™ ax ~ l 
is described by the representation T to** -1 , which according to the above coincides with 
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Ts m ax on 7i^ ax ~ 1,Smax which is therefore an invariant subspace. Therefore T s ™ ax ~ l is not 
irreducible therefore it is fully reducible (because it is a finite dimensional representation 
of a finite group). Therefore the orthogonal complementer of 7if™ aa: ~ 1 '' Sma:E which is the 
f-lS maa .-i,s max -i j g a j go an i nvar j an t subspace. The restriction of T 5 ™^ 1 to this subspace we 
denote by T Smax -\ := T s ™ a *~ x \ s^ ax -i,s max -i- Continuing this process on this line, acting 

f*-s 

by S~ again and again, we get: T s * = ®g=, Sie , T s . 

In other words this means that the functions <^f z ' , which describe the spin configuration 
with S and S z values of the total spin and its z component for fixed hole configuration, form 
a space of a representation of the symmetric group. This above defined representation 
depending only on the value of S was denoted by T s , and the effect of the Hamiltonian on 
the spin configuration can be described by this representation: 

Hg\ V T' S ) = T S [C(TZ^^)}\ V T' S ) (28) 

The matrices of the representations Ts can be given by algebraic methods, see e.g. 
Hamermesh 1962. These methods are usable only in case of few particles. In this case we 
can use the expression Eq.fllOD to get exact diagonalization results. 
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FIG. 1. The a: sqare of the total spin and b: specific heat as a function of the reduced 



temperature in the case of 5 x 5 lattice with two holes. The calculation based on the Eqs.17-18 
in which the coefficients N^ c \l) are taken into account exactly up to I = 14 and determined by 
Monte Carlo simulation up to I = 50. 
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